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We present a method to calculate integrals over monomials of matrix elements 
with invariant measures in terms of Wick contractions. The method gives 
CNj ' exact results for monomials of low order. For higher-order monomials, it 

leads to an error of order 1/iV" where A'" is the dimension of the matrix and 
^ . where a is independent of the degree of the monomial. We give a lower bound 

^] , on the integer a and show how a can be increased systematically. The method 

\ is particularly suited for symbolic computer calculation. Explicit results are 

' given for 0(A^),U(A^) and for the circular orthogonal ensemble. 
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I. INTRODUCTION 



The calculation of group integrals of monomials of matrix elements for compact Lie 
groups has a long tradition going back to Ullah and Porter Their results were later 
extended to the space of symmetric unitary matrices endowed with Dyson's invariant 
measure P]. The problem was nearly dormant for some years but was recently solved 
completely for the orthogonal group by recursion 0]. Unfortunately, it seems that there is 
^ I no easy generalization of this method to other groups. The method proposed in Ref. 0, 
I on the other hand, is quite general but soon becomes rather cumbersome. Moreover, that 
method is not suited for computer-supported analytical work. 

Aside from their immanent group-theoretical significance, group integrals over mono- 
mials of matrix elements for compact Lie groups are important in applications of random- 
matrix theory. This field has seen an explosive growth over the last decade and has become 
an important tool in a great variety of fields in physics, chemistry and related areas 0. This 
fact lends further urgency to the evaluation of the above-mentioned group integrals. 

In this paper, we present and analyse a novel approach to the problem. We encountered 
the problem in the study of a random-matrix model for a class of chaotic systems (the 
semi-separable systems) p. We found that we could use the invariants of the orthogonal 
group to construct a weight function. With the help of this function, it was possible to 
evaluate the group integrals in question by simple Wick contractions. We implemented the 
scheme on the computer and found that beyond the expected exact results for low-order 
monomials, monomials of higher order were also calculated correctly up to and including 
the sub-leading order in where N is the dimension of the matrices under consideration. 
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We conjectured that this statement holds for monomials of any order. It is the purpose of 
the present paper to extend and prove the conjecture and to explore the scope of its validity 
beyond the orthogonal group. 

The group integrals extend over a compact matrix space with a measure which is uniquely 
determined by the underlying symmetry group. In order to reduce the computation to Wick 
contractions, we consider an extended matrix space where all matrix elements are indepen- 
dent Gaussian variables. In this space, all integrations trivially reduce to Wick contractions 
and can easily be implemented in many programming languages. The constraints due to the 
group structure are then introduced in an approximate fashion through a weight function w 
appearing as factor in the integrand. This function is chosen in such a way that the integrals 
yield the exact values for the lowest-order invariants of the group. It turns out that w is 
not always positive and, thus, not a measure. This, however, is not a significant obstacle. 

Candidates for our spaces are the orthogonal group O(A^), the unitary group U(A^) = 
CUE(iV), also known as the circular unitary ensemble, the circular orthogonal ensemble 
of symmetric unitary matrices COE(A^), and the circular symplectic ensemble which is 
isomorphic to the unitary symplectic group CSE(A^) = USP(2A^). We note that COE(A^) is 
not a group — for all other cases we can use the Haar measure, while in this case we have 
to use Dyson's invariant measure. 

We first present our arguments for the case of the orthogonal group. In Section y we 
construct the weight function w from the invariants of the orthogonal group. We show that 
the defining equations for w always have a unique solution. We give explicit expressions for 
w in the simplest cases. We show that low-order monomials are calculated exactly using 
Wick contraction. In Section |IT|, we show that monomials of higher order are evaluated 
correctly by Wick contraction, up to an error of order A^~". We establish a lower bound for 
the exponent a. In Section ^ we extend our arguments to other matrix spaces. We give 
explicit expressions for the weight functions for U(A^) and for COE(A^). The more involved 
and less important case of CSE is only touched upon. 

II. THE WEIGHT FUNCTION FOR 0{N) 

As explained in Section |, we start with a space of real matrices M. The elements are 
taken as independent Gaussian-distributed variables with zero mean and identical variances. 
In other words, our measure dfig for integration is the product of the differentials of all matrix 
elements times N" exp{—N trace(MM-^)} where A/" is a normalization factor. We recall that 
N is the dimension of the matrices M. We are interested in values of ^ 1. The measure 
is invariant under right or left multiplication of M with any orthogonal matrix. To restrict 
the integration to the orthogonal group, we could think of multiplying dfj,g with a set of delta 
functions expressing all the constraints due to orthonormality. This is clearly impractical. 
Instead, we modify the measure by multiplying dfig with a weight function w^. This function 
is chosen in such a way that all orthogonal invariants up to and including order 2k are exactly 
reproduced when we use as a weight function under the integrals over M. We note in 
passing that our present notation differs from that of Ref. p. Our index k equals two times 
the index k used there. 

To determine w^, we consider all invariants lj{0) up to order 2k in the matrix elements 
O of the orthogonal group. Here, j is a running index. We recall that all such invariants are 
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even in the O's. In every invariant lj{0), we replace O by M to obtain Ij{M). We write w^, 
as a linear combination of all the Jj(M)'s up to order 2k in M. The coefficients of the linear 
combination are determined by the requirement that the average of every Ij{M), calculated 
by integration over dfig with weight function w^, yields the same result as integration of 
lj{0) over the orthogonal group. We recall that with n a positive integer, the invariants of 
the orthogonal group are given by expressions of the form tr{(00^)"}, or by products of 
such expressions. We accordingly write the quantities Ij{M) in the form 

jW(M) = n tr{(MM^)'='} . (1) 

Here 2k denotes the degree of / in M, and k = {ki, k2, . . .) is a partition of k into positive 
integers k^ > 1 with ki+k2+. . . = k. Without loss of generality we require that ki > k2 > . . .. 
The weight function is now written as 

w.iM)=ai;'^ + ±j:ai^hl'\M). (2) 

fc=l k 

The sum on the right-hand side of Eq. (j^) extends over a complete set of linearly independent 
invariants up to order 2k in M. 

(k) 

We determine the coefficients from the conditions of orthonormality. More precisely, 
we require that the relations 

I dfigW^iM) = 1, 

J diigW^{M){MM^)i,,, = 6i,j, , 
J dfigW^{M){MM^)i^,^{MM'')i,^, = 6i,,,,6i,j, , 

J dfigW^{M){MM^)i,^^ {MM^)i,j, . . . {MM^)i^,^ = 6i,j,6i,j, x . . . x 5^^,,^ (3) 

be fulfilled identically. Relations of the form (^) hold for any value of k for the orthogonal 
group but must be imposed for the integration over the matrices M. 

Eqs. determine the coefficients a|f^ uniquely. To show this, we take traces over 
these equations in such a way that the integrals on the left-hand sides take the form 
/ djjLgW;^{M)l'^\M). The resulting set of equations has the form 

J dfigW,{M)li'\M) = bP (4) 

where the coefficients B^^ are given by powers of A^, with N the dimension of the matrices M. 
Recalling Eq. (|), we see that Eqs. (|) constitute a set of linear equations for the coefficients 
There are obviously as many equations as there are coefficients We conclude that 
Eqs. d^) possess a unique solution unless the determinant of the matrix C with elements 
^kiifa'^ = / dfiglj^_^^ ij^^^ vanishes. But if det(C^^^J^^^'') = 0, there exists a nontrivial solution 
b^^^ of the homogeneous equation I]fc2k2 ^kikt^'^^kt^'' — 0- The existence of this solution 
imphes that we also have J2kik2 I^kik2 ^ki^ ^^1^2^ ^^2^ = 0- Recalling the definition of the 
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matrix C, we observe that the last relation can be written as / dfig\ J^ki^i -^ki^''^£^^P = 0- But 
the integrand in the last expression is positive semidefinite and does not vanish identically. 
Therefore, it is not possible that det{Cl^^^^^) vanishes, and the solution of Eqs. (H) exists 
and is unique. This solution also solves Eqs. (|^). To see this, let us assume the contrary and 
focus attention on the second of Eqs. (^. (The argument is easily extended to the entire 
set of Eqs. (|^)). Inserting the solution of Eqs. (Q) into the left-hand side of that equation 
yields on the right-hand side the terms + Aj^j^ where the matrix A is both traceless 
and invariant under every orthogonal transformation. This implies A = 0, in contradiction 
to the assumption that we did not find a solution of the second of Eqs. (^). 

Eqs. (§) imply that the integrals over all polynomials of degree n < 2k in M are equal to 
the corresponding expressions for O(A^). To see this, it suffices to consider the integral over 
an arbitrary monomial of degree n. It is obvious that the integral vanishes unless n is even, 
n = 2k. We write the monomial as A4^"'^ = Mi^j^Mi^j^ . . . Mi^j^. The integral over Ai^"'^ is 
obviously invariant under right or left multiplication with any orthogonal transformation. 
Therefore, the integral over Al^") must be a linear combination of invariants multiplied by a 
suitable set of Kronecker deltas in the indices ii, . . . , i„ and ji, . . . , j„. By construction the 
invariants have the same values as in O(A^). 

Inspection shows that the weight function wq = 1 fulfills the second of Eqs. (^) automat- 
ically. Thus, Wi = Wo and, therefore, ag^^ = 1, a^i^ = 0. The first nontrivial condition is, 
therefore, the one appearing in line 3 of Eqs. (|^). This condition (and all that follow below 
it) is violated by Wq. We now give the explicit results for the first few weight functions w^- 
These were obtained with the help of the Mathematica program. For k = 2, we find 
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" ~ 4(-l + iV)(2 + iV) ■ 



For K = 3, we have 
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5N^ 



A{-l+N){2 + N) 

+ N){2 + N) 
ATS 



3(-2 + N){-1 + iV)(2 + A^)(4 + N) 
~ (-2 + iV)(-l +N){2 + N){4: + N) 
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(3) 2iV3 



3{-2 + N){-l + N){2 + N){'i + N) ' ^ ' 
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(4) 
021 



"111 



(4) 
O2II 



For K — A, we have 

-60 iV^ + iV^ 



16 (-1 + iV) (2 + AT) 

56 iV2 + 2 + 
16 (-1 + N) {2 + N) 
7N^ 



3 (-2 + N) (-1 + iV) (2 + iV) (4 + N) 

-AS N"^ - 2 - 
8 (-2 + A^) (-1 +N) (2 + N) (4 + A^) 

88 A^=^ + 6 AT^ + 3 A^5 
24 (-2 + A^) (-1 + A^) (2 + A^) (4 + A^) 

N'^ (6 + 5N) 



"31 — 



"22 



8 (-3 + N) (-2 + N) (-1 + N) {1 + N) (2 + AT) (4 + AT) (6 + AT) 

A^6 (6 + 5 A^) 

(-3 + A^) (-2 + AT) (-1 + A^) (1 + N) (2 + AT) (4 + A^) (6 + A^) 

AT^ (18 + 5Ar + Ar2) 
32 (-3 + N) (-2 + N) (-1 + AT) (1 + A^) (2 + A^) (4 + N) (6 + N) 

{72 + 78N + 5N'^ + N^) 
~ 16 (-3 + A^) (-2 + N) (-1 + AT) (1 + A^) (2 + A^) (4 + A^) (6 + A^) 



(4) A^^ (72 + 78A^ + 5Ar^ + AT^) 

"^^^^ ~ 32 (-3 + A^) (-2 + A^) (-1 + A^) (1 + A^) (2 + A^) (4 + A^) (6 + AT) . ^ ' 

We note that with increasing k, the expressions become rather involved. Moreover, the 
coefficients a^^ with the same lower indices k change with k. 



III. MONOMIALS OF HIGHER ORDER FOR 0{N): THE 1/N EXPANSION 

We have seen that integrals over all polynomials of degree n <2k have the same values as 
for O(A^). What about polynomials of higher order? Again, it suffices to consider monomials 
of even degree 2k with k > k. We show that the integral over A^*^'^) coincides with 
the result for O(A^) up to terms of order A'^"" where the integer exponent a is positive and 
independent of k. More precisely, we show that for k > k, we have 

2k 2k 

N'^ [ diigW,{M) n M,,,-, ^N^ f dho^N) n Oi^j^ + ^ (1/^") where a > [k/2] + 1 . (8) 

•' V=l ^ I/=l 
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Here [/t/2] indicates the integer part of k/2, and (i/io(Af) denotes the Haar measure for 
integration over 0{N). We note that the factors A^^ in front of the integrals normahse the 
A^-dependence so that these terms are (at most) of order 1. Another equivalent form of 
Eq. (H) is obtained by summing over pairs of indices ji = j^, js = ji, ■ ■ ■■ This removes the 
factors A^'^ and yields 

k k 

f dfigW^iM) n {MM^),^i^ = J dhoiN) n (OO^^i^ + O (l/AT") where a > [/t/2] + 1 . 

(9) 

The equivalence of Eq. (P) and Eq. @ follows from the fact the matrix C discussed above, 
if defined with respect to properly scaled monomials, does not depend on A^. The remainder 
of this Section is devoted to proving Eq. (^. 

It is useful to introduce a few auxiliary concepts. We consider Gaussian integrals over 
monomials of M without the weight function w^. We write for brevity 

k 

J df,, n (MM^j^ = {{MM^f), (10) 

where the index g indicates the purely Gaussian integration. To define the completely 
correlated part of this expression, we consider first the case k = 2. We use Wick contraction 
and have 

{{MM^Y), = {{MM''))g{{MM^))g + ((MM^)2),, . (11) 

The last term on the right-hand side of Eq. (0) is the completely correlated term. For 
the general case of arbitrary order 2k, we define the correlated part {{MM'^)^)gc as that 
contribution to {{MM'^)^)g which cannot be written in the form of products of two or more 
factors, each of which is a complete Wick contraction of powers of MM'^. It is easy to see 
that 

{{MM^)'')g, = 0(l/N''~^) . (12) 

The linear increase with k in inverse powers of A^ in Eq. (|T2D is due to the fact that every 
Wick contraction which connects two M's appearing in different factors MM'^ suppresses 
one summation index. Therefore, the correlated part {{MM'^)^)gc contributes the highest- 
order terms in 1/A^ to {{MM'^f )g. 

We now consider integrals involving the weight function and use the same notation, 

k 

J dfigw^ n iMM^),^i^ = {w,{MM'^f)g . (13) 

Again using Wick contraction, we define the correlated part {w^{MM'^Y) g^ of this expression 
as that part which cannot be written as the product of two or more factors, each of which 
is a complete Wick contraction of powers of MM^ and/or w^. 

We proceed to show that in the equations relating the integral {w,i{MM'^)^)g to the 
integral over the Haar measure, the leading correction term (lowest order in 1/A^) which 
does not cancel is given by 
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{w^{MM^f)g, = O (l/Ar[('=+i)/2]) for l<k<K . (14) 

This relation is based upon the assumption that there is no accidental cancellation among 
the terms contributing to lowest order in Therefore, [{k + l)]/2 actually constitutes a 
lower bound on the exponent of 

To prove the relation (|T3p, we rewrite the defining equations for w^, Eqs. (^, as follows. 
We consider the expression {wi^{MM'^)^)g with k integer and k < k. We decompose this 
expression into correlated contributions. These originate from all partitions k = {ki, k2, . . .) 
of k with ki > fcj+i and J2i h = k. We denote by io the smallest index for which all kt with 
i > io are equal to one. Then, we have 

k i k^l'= i jj^i 

In the first term on the right-hand side, we have used that {wi^)g = 1. In the second term, 
we have used that {Wfi{MM'^))gc = 0. Trivially, the second sum on the right-hand side of 
Eq. (|15D does not extend over the partition 1^ = (1, 1,1,.. .) {k terms unity). According to 
Eqs. (^), the expression in Eq. equals (((MM^))^)'^. This equals the contribution from 
the first sum on the right-hand side for the partition l'^. The remaining terms must vanish, 

io 

E IlmM^)'')9c + E T.{MMM^r^)g^l[{{MM^r^)g^ = . (16) 

Eq. ([T6|) must hold for all values of k with k < n. To proceed, we observe that the partitions 
of k can be grouped into classes as follows: Partitions within the same class carry the 
same number p of /cj's that have value unity. The classes are labeled by p, namely Cp. For 
instance, for k = 6, class C2 contains the partitions (4, 1, 1) and (2, 2, 1, 1). In Eq. (|16D , we 
order the sum over k by grouping together all partitions which belong to the same class. 
We show presently that each such contribution must vanish separately. Then, we have for 
every p = 0,1, . . . , k — 2 that 

io 

E ]lmM^)'')9c+ E T.{MMM^)>'')g^l[{{MM^p)g^ = . (17) 

k7^1,k6Cp keCp i jf^i 

Eq. ([TtD follows directly from the facts that Eq. ([TBI) holds for all k < k, and that the 
contributions from class Cp to a partition of k are the same as the contributions of class Cq 
to a partition of k—p, except for a string of Kronecker delta's due to the factors {{{MM'^)) gY . 

We are now in the position to prove the relation (|T^). We observe that in Eqs. (|T3), the 
term {w^{MM^Y)gc appears only in the class Cq. Therefore, we have 

{W^{MM^f)g^ = - E Y[mM^t)g^ 

kgCo i 

- E Y.(MMM^t)g^\[{{MM^p)g^. (18) 

kGCo,k7^(fc) i jjti 



Using complete induction, i.e., assuming that the relation ([T^) holds for all values of k' with 
k' < k, and using Eq. (|12]), we conclude from Eq. (|18D that the relation ([I^) also holds for 
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k' + 1 = k. The terms of lowest order in originate from partitions which have either 
the form (2, 2, 2, . . .) (for even k) or (3, 2, 2, 2, . . .) (for odd k). Again, we cannot rule out 
the occurrence of accidental cancellations which would increase the power [{k + l)]/2 in the 
relation (||). 



Having established the relation (p^, we turn to the center piece of this Section, Eq. (^. 
We first consider the case k = k + 1 and decompose the integral into correlated terms with 
contributions from all partitions of k + 1. For these partitions, we write K + 1 = [ki, /c2, . . .) 
with ki + k2 + . . . = K + 1 and ki > k2 > ■ ■ ■■ We use the notation introduced above. Then, 

P (K+l)eCp i 

+E E E(^^(^^^)'0.cn((^^^)'^).c- (19) 

P (K+l)eCp i jjti 

Eqs. ([T7| ) imply that all terms with p ^ and p ^ k + 1 vanish, and we are left with 

{w,{MM^r+')g = (((MM^)),)-+i + E I[(iMM'')'')gc 

(K+l)eCo i 

+ E J:{Mmm^)'^),^yI{{mm^P),^. (20) 

(K+l)eCo * jT^i 

We use the same argument as in the previous paragraph and the result (|T^ and Eq. (|r^) . We 
conclude that the contributions of lowest order in result from the partitions (2, 2, 2, . . .) 
or (3,2,2,2,...), respectively, and arrive at Eq. (^. The exponent a has the value a > 
[k/2] + 1. In writing the inequality sign, we again allow for the possibility that an accidental 
cancellation of contributions from these partitions occurs. 

We finally have to consider the case where k — k = n > 1. Decomposing the expression 
{Wfi{MM'^Y^^)g as in Eq. (|19]), we easily see that the terms of lowest nonvanishing order in 
stem from the partitions in class C„_i. But these give exactly the same contributions 
as those for = k + 1 that were estimated in the last paragraph. This completes the proof 
ofEq. ®. 



IV. THE UNITARY GROUP AND OTHER UNITARY ENSEMBLES 

In this Section, we primarily address integrals over momomials of unitary matrices U with 
respect to the Haar measure of the unitary group. We start with Gaussian integrals over 
complex matrices M. The real and imaginary parts of the matrix elements are independent 
and Gaussian-distributed. The integrals are again worked out using Wick contractions. For 
the weight function w^, we write in analogy to Eq. (Q) 



<{M) = h\^^ + Y.Y.^'C^t\M). (21) 

A;=l k 



The invariants are defined as in Eq. (|I]) with MM'^ replaced by MM'^ where ^ stands for 
Hermitean conjugation. The upper index u refers to the unitary case. The values of the 
coefficients 6£ are, of course, not the same as for the orthogonal case. The arguments for the 
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existence and uniqueness of the solutions carry over without change, and again computer 
programs are available to perform the contractions and calculate the coefficients b^. For 



and we find 



and 



bt 



b'n 



bl 



bt. 



bt. 



-'22 



bin 



& = 1 - — 



bl 
bl 



bl. 



1 

N 



Ar2 
2" 



Ar3 



2(A^- l)(iV+ 1) 
2(A^- 1)(A^ + 1) ' 



(22) 



24 - 46 + 3 
24 

;-i2 + iv2)) 



A^=^ (-30 + N^) 



4 



-^1111 



-1 + N) (l + N) 
(28 + N^) 



-1+N) (l + N) 



-2 + N) (-1 + iV) (l + N) {2 + N) 
-iV^ (24 + N^) 



-2 + N) (-1 + AT) (1 + N) (2 + N) 
(44 + 3 A^2) 



~2 + N) i-l + N) (l + N) (2 



N) 



-3 + N) {-2 + N) 



-1 + A^) (1 + N) (2 + A^) (3 + N) 
5N^ 



3 + N) (-2 + A^) (-1 + N) {1 + N) (2 + A^) (3 + N) 
N^ (6 + N'^) 



-3 + A^) (-2 + A^) (-1 + A^) (1 + N) (2 + A^) (3 + A^) 
-N^ (36 + N^) 



-3 + A^) (-2 + AT) (-1 + AT) (1 + N) (2 + A^) (3 + A^) 
A^^ (36 + N^) 



[-3 + N) (-2 + A^) (-1 + N) {1 + N) (2 + A^) (3 + A^) ' ^^^^ 

The arguments determining the leading contribution to the 1/A^ expansion are the same 
ones as for the orthogonal group. Thus, we have to increase k by two to improve the error 
by one order in 1/A^ in the calculation of monomials of high order. This is the reason for our 
not giving but only which yields correct values for the integrals up to order l/N"^. 

For COE and CSE, the situation is slightly more complicated. The constraints on the 
matrices are not expressible in a simple way in terms of products as done in Eqs. (|^). It 
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seems, therefore, most convenient to limit the space of independent matrix elements from 
the outset. In the case of the COE this is fairly simple: The symmetry reduces the number 
of independent complex matrix elements to N{N + l)/2. Therefore, we consider Gaussian 
averages in the space of complex symmetric matrices S = S"^. We accordingly have a 
contraction rule with two terms, 

J dfig{S)S*jSki = {SikSji + SiiSjk) ■ (24) 



The invariants are now defined in terms of SS*. Introducing the rule (p^) into our program 
we can again calculate w^. As an example we the coefficients for W2, 



4 = 1- 



N{N + 1) 



-11 



N + 1 
2 

{N + lf 
4:N{N + 3) 

AN{N + 3) ■ 



(25) 



In addition, the following subtle point must be considered. In Section ^ we have used 
the invariance of the Haar measure to justify that contraction with gives exact results 
for all polynomials up to order 2k in the matrix elements. In the present case we have no 
invariance group and by consequence no Haar measure. On the other hand Dyson's measure 
with respect to which we integrate, is also totally defined by an invariance group albeit a 
smaller one than that of U(A^). The important point is that again the measure is uniquely 
defined by a linear group of transformations. The orthogonality conditions resulting from 
the unitarity of the matrices are the same and symmetry is taken into account explicitly in 
the contractions. Therefore all arguments again go through and indeed inspection of the 
results obtained by our code with those obtained in Ref. shows agreement. 

The case of CSE is simpler because it involves an invariance group, but more compli- 
cated because the matrices are symplectic. Two ways seem open to address this case. We 
might include the symplectic property from the outset in the contraction rules, or we might 
introduce this property as a constraint in the expression for the weight function. While both 
ways seem possible it is not clear which one is easier to follow. In view of the fact that CSE 
is of minor importance for practical applications, we have left this problem open. It is clear, 
however, that it can be tackled along the same lines. 



V. CONCLUSIONS 

We have presented a systematic way to calculate integrals over monomials of matrix 
elements for compact matrix groups and for other matrix ensembles whose measure is de- 
fined uniquely by an invariance group, such as the circular orthogonal ensemble of unitary 
symmetric matrices. This method gives exact results for monomials of low order. For higher- 
order monomials, it leads to an error of order which is independent of the degree of 
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the monomial. Wc have given a lower bound on the integer a. and we have shown how a can 
be increased systematically. The method is particularly suited for symbolic computer calcu- 
lation. Codes are available for 0(iV), U(A^) as well as for the circular orthogonal ensemble 
in Mathematica and in C from one of the authors (T.P.). 
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